Andreev spectra and subgap bound states in multiband superconductors 
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A theory of Andreev conductance is formulated for junctions involving normal metals (N) and multiband 
superconductors (S) and applied to the case of superconductors with nodeless extended s±-wave order param- 
eter symmetry, as possibly realized in the recently discovered ferropnictides. We find qualitative differences 
from tunneling into s-wave or d-wave superconductors that may help to identify such a state. First, interband 
interference leads to a suppression of Andreev reflection in the case of a highly transparent N/S interface and to 
a current deficit in the tunneling regime. Second, surface bound states may appear, both at zero and at non-zero 
energies. These effects do not occur in multiband superconductors without interband sign reversal. 

PACS numbers: 74.20Rp,74.50.+r,74.70.Dd 
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The recent discovery of high-T c superconductivity in fer- 
ropnictides has been a major event in solid state physics. The 
first theoretically proposed pairing symmetry for this com- 
pound, s-wave with sign-reversing order parameter (s±) has 
been followed by a number of theoretical papers substantiat- 
ing this proposal at various degrees of sophistication and ex- 
ploring ramifications of the proposed state 01 . Within this 
proposal, two sets of Fermi surfaces are distinguished: the 
hole Fermi surfaces around T and the electron Fermi surfaces 
around M. The it phase difference between the superconduct- 
ing condensates is thought to be induced by spin-fluctuations. 
Experimental evidence has been favorable to the s± model so 
far, but is still ambiguous. 

Andreev spectroscopy is a strong experimental probe of the 
superconducting order parameter. But in the case of the fer- 
ropnictides, both nodeless as well as nodal superconductiv- 
ity have been inferred from the absence JUl or presence J]!] 
of zero-bias conductance peaks. Also, point-contact spec- 
troscopy has both provided evidence for a single gap as well 
as multiple gaps J2, 01. One of many questions to be asked 
in this connection is how possible interference between the 
two bands in the ferropnictides may affect the Andreev con- 
ductance spectra. Can the interference phenomenon be used 
to distinguish the s± state from other scenarios? To address 
these questions, a generalized theory of Andreev conductance 
is needed, relevant also to other multiband systems. 

Surface phenomena in s± superconductors have attracted 
considerable recent attention El H S Si ■ Certain limiting 
cases were considered, but no general calculation of Andreev 
and tunneling current, properly accounting for the effect of 
interference between the two relevant bands, have been pub- 
lished so far. This we provide in the present Letter. The for- 
mation of bound states at a free surface of an s± superconduc- 
tor, at an S±/N/S, and at an N/S/S± junction was found in 
Refs. |4J, |5| and la, respectively. However, the conditions for 
such a bound state and its effect on Andreev and tunneling 
conductance were not addressed in these papers. Ref.0 found 



an enhancement of the density of states at zero energy in a 
thin N layer on the top of an s± superconductor, but it is un- 
clear whether this numerical result is related to bound states 
or not, since finite energy rather than zero-energy bound states 
were predicted in J3, H 01 ■ Finally, Ref . [H considered the 
same problem as ours, the conductance spectra in an N/S± 
junction, but did not find any new effects compared to regu- 
lar Andreev reflection. This may be related to not properly 
accounting for any interference effect. At the moment, a gen- 
eral analytical unrestricted treatment of an interface in an ar- 
bitrary s± superconductor seems necessary to clarify the ex- 
isting confusion. 

We have studied Andreev conductance in an N/S± junc- 
tion by including in the classical "BTK" model fgl the inter- 
ference between the excitations in a two-band superconduc- 
tor at arbitrary interface transparency. Applying our extended 
BTK model to the s± scenario we have found qualitatively 
new effects. For a fully transparent interface, the destructive 
interband interference leads to a strong suppression of An- 
dreev reflection, in striking contrast to the conventional case. 
In the tunneling regime, two new effects are found: (a) a 
current deficit at high bias voltage, which is also due to de- 
structive interband interference, and (b) Andreev bound states 
(ABS) similar to those responsible for the zero-bias anomaly 
(ZBA), a well known fingerprint of d-wave superconductors 
111 Oil . However, instead of a ZBA, a peak of similar origin 
may appear (depending on the parameters) at a finite energy, 
that can easily be mistaken for an extra gap. 

A ballistic Andreev contact can be modeled by a one- 
dimensional conductor, whose right half (x > 0) is a two- 
band metal (two different states at the Fermi level, one with 
the wave vector p and the other with q), and the left half is 
a simple metal. The wave function at the energy Ep in the 
left half is ^ = ipk(x) + bip-k{x), where the first term 
is the incident Bloch wave and the second term the reflected 
Bloch wave. The wave function on the right-hand side is 
^r(x) = c[cf)p(x) + ao<p q (x)]. Here p and q are the Fermi 
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vectors for the two bands, cf> is the Bloch function in the two- 
band metal, and the mixing coefficient ccq defines the ratio 
of probability amplitudes for an electron crossing the inter- 
face from the left to tunnel into the first or second band on 
the right. Similar problems have been studied in the theory of 
the tunneling magnetoresistance, where the leads are usually 
multiband d-metals. 

The main conceptual pitfall here is that the standard ap- 
proaches to tunneling assume the wave functions to be plane 
waves. However, there cannot be two different plane waves 
propagating in the same direction with the same energy. This 
may only occur when the wave functions are Bloch waves 
- which they are in reality. It has been realized in the last 
decade that there is not a single factor that would define the 
relative tunneling probability of the two Bloch waves. If the 
wave vector parallel to the interface is not conserved (it usu- 
ally is not, except perhaps for epitaxially grown contacts), one 
factor is the number of tunneling channels in each band, pro- 
portional to (Nv±), where N is the density of states and v± 
is the Fermi velocity component normal to the interface. Even 
more important is the character of these Bloch wave functions. 
E.g., states of different parity on the right and on the left sides 
of the interface hardly overlap, so that even a weak interface 
barrier will strongly suppress tunneling from particular bands. 
With this in mind, we keep ccq arbitrary and present the results 
for different cases. One implication is that the observable tun- 
neling spectra may actually change drastically from contact to 
contact, as the interface properties change. Indeed, there are 
indications that this may be the case 111 ill . 



ary conditions at the interface (x = 0) 



At a normal metal (N) - superconductor (S) contact, in the 
case of the two-gap model with unequal s-wave gaps, one can 
write 



*(0) = * s (0) = 4-^(0), 



(3) 



£^°)-£^°)=^' (4) 

where H is the strength of the (specular) barrier. 

The boundary conditions on the wave function derivatives 
are usually expressed in terms of Fermi velocities. However, 
a closer look reveals that this is actually incorrect for Bloch 
waves. Therefore, in the following, we introduce an "inter- 
face velocity". For a given Bloch function, say, ipk(x) = 
Y^q Ac,k exp[i(fc + G)x], it is defined as 



Vk 



ih~ 1 dipk(x) 
mijjk(x) dx 



\x=0- 



(5) 



The so defined Vk is real and has the same symmetry prop- 
erties as the Fermi velocity (this can be shown by expanding 
the wave functions in terms of the plane waves), but it coin- 
cides with the actual group velocity only for free electrons. 
For general Bloch waves it is different, and even dependent 
on the position of the interface plane in the crystal. We leave 
the interesting and important issue of the relationship between 
the interface velocity and the group velocity Il2ll for a further 
study, and proceed with the problem at hand. 

Now, introducing the barrier strength Z = H/Hvn, where 
vn is the velocity on the N side, defined according to Eq. ©, 
and using the boundary conditions Eqs. ([3] IU, we find the 
general solution for a, b, c, and d. It depends on Z and on 
the ratios of the interface velocities. To keep the expressions 
compact, we present them below for the case of equal inter- 
face velocities on the N side and in both bands on the S side. 

For the s± model where tpi — ip2 = tt this gives 
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Here tp 12 are the phases of the gaps Ai ; 2 in both bands, 
u and v are the standard Bogoliubov coefficients u\ 2 = 

i + ^E 2 - A 2 2 /2E, v 2 2 = i - Je~ 2 - Aj 2 /2E. In the 
case of the s± gap model with unequal s-wave gaps of oppo- 
site sign we have ipi — tp2 = 7r, while in the standard two- 
band model with the gaps of the same sign we have ip\ = (f2- 
The amplitudes a, b describe Andreev and normal reflection, 
and the amplitudes c, d describe transmission without branch 
crossing and with branch crossing, respectively. 

The total wave function must satisfy the following bound- 



7a = uivi — a{u\V2 + U2V1) + a 2 U2V2, 
76 = (Z 2 + iZ) [v{ -u\ + ol 2 (u\ - v\ 
7c = (1 — iZ)[u\ — au2)5, 
jd = iZ{v\ — av2)S, 
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where 7 = (1 + Z 2 )(uj - o?u\) - Z 2 {vj - a 2 v 2 ), S = 
ipk(0) / <t> P (0) and a = a O (j) q (0) / 4> p (0) . Note that for plane 
waves <5=1 and a = 0. 

For the ,s ++ model with ipi = ip2 we obtain 

7a = U\V\ + a(uiv 2 + u 2 vi) + a 2 u 2 V2 , 

76 = (Z 2 + iZ) [(ui + av 2 ) 2 - (ui + a« 2 ) ] , 

7c = (1 — iZ)(m + CW2) 5, 

yd = iZ{v\ + av 2 ) S, (7) 

with 7 = (1 + Z 2 )(u x + au 2 ) 2 - Z 2 (v 1 + av 2 ) 2 - 

In a single band case (a = 0) and for plane waves the stan- 
dard BTK results are recovered. Below we shall discuss new 
effects arising in the s± model. First, we consider a trans- 
parent interface, Z = 0. In the s± case, b = d = 0, 
a = (vi — aV2)/(ui +01112), c = l/(i*i + 01U2). At E = we 
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get a = (VA^— av / A^)/(v / A ; j~+ ay/A^) < 1, i.e. Andreev 
reflection is suppressed. On the other hand, in the s ++ case 
b = d = 0, a = (vi + a,V2)/(ui + ctU2), c = l/(«i + otu^) 
resulting in a = 1 at zero energy, as expected in the stan- 
dard BTK model. This effect is due to the destructive interfer- 
ence between the transmitted waves in the s+superconductor, 
which was missing in the previous works JJlHlaSISl- 
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, then the probability current Jp is given by 



J P = - [Im(/*V/) -Im [g* Vg)] 
m 



(8) 



properly taking electron and hole contributions into account. 
Using Eq. (f2]) for * at the superconducting side of the in- 
terface, Jp = (C + D) Jk, where Jk = vn ^-(O)! 2 is the 
probability current of a normal electron in the state ijjk, and 
the transmission probabilities C and D depend on the veloci- 
ties in the two bands. For equal band velocities they are given 
by 



C = \c/5\ 2 [wi + a 2 w 2 + 2aRe (u lU * ± v^)] , (9) 
D = \d/5\ 2 [wi + a 2 w 2 + 2aRe ( Ul u* 2 ± viv%)] , (10) 

s ++ models respectively, where Wi 2 = 



for the s± and 

\U1,2\ 2 ~ |^1,2 | 2 - 

At the normal side of the interface the probability current 
is (1 - A - B) J k , where A = \a\ 2 and B = \b\ 2 are An- 
dreev and normal reflection probabilities. From Eqs. (6,7) 
and (9,10), it can be verified that A + B + C + D = l. Thus 
we have proven that the probability current is conserved. The 
electric current / across the contact is given by the standard 
expression |@] 
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[/o (E - eV) - f (E)] [1 + A-B] dE, 



(11) 

where /o is the Fermi function, Rjy is the normal state inter- 
face resistance and V the voltage bias across the interface. 
Below, we present the results of calculations of the angle- 
resolved conductance dl/dV in the regime of a fully trans- 
parent interface Z = and in the tunneling regime Z ^> 1. 

The zero-temperature conductance at Z = is shown in 
Fig. 1. In the s++ case, there is a standard enhancement of 
conductance at low bias eV < Ai due to Andreev reflection, 
followed by a decrease of conductance and a weaker feature 
at eV = A.2- At the same time, as discussed above, a strik- 
ing suppression of the zero-bias conductance occurs in the s± 
case. The strongest suppression occurs at a — y A1/A2. 

Figure 2 shows the zero-temperature conductance for large 
Z in the s± case. Sharp conductance peaks appear at certain 
values of a. These peaks have a clear interpretation as An- 
dreev bound states. Indeed, for large Z, a bound state exists if 
7 = 0, that is, if u\ — v\ = c?{u^ — 
bound state is 
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FIG. 1 : (Color online) Conductance in the case of fully transparent 
interface, Z = 0: comparison of the s± and s++ models. 
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FIG. 2: (Color online) Conductance in the low transparency regime, 
Z — 10, in the s± model. 



If Ai = A2 this gives the trivial Eb = A solution, that is, no 
subgap bound states. If a = 0, similarly, Ep = Ai. However, 
when < a 2 < A1/A2 bound state solutions exist (see Fig. 
3), most notably a zero-bound state Ep — if a 2 — A1/A2. 
Note the bound states for a = 0.5 and 0.7 in Fig. 2. 

Further, it is also seen from Fig. 2 (e.g., for a = 0.9) 
that there is a current deficit at high bias. This feature is 
due to a destructive interband interference and is in contrast 
with the properties of N/S junctions known so far, irrespec- 
tive of whether S is s- or a d-wave. The only known case is 
a double-barrier junction, where current deficit occurs due to 
non-equilibrium quasiparticle distribution in the N layer 1 1311 . 

For comparison, we also present the results for the s ++ 
model in Fig. 4, where bound states are absent, as expected. 
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we have concentrated on the discussion of the angle-resolved 
conductance. The total conductance depends on the orienta- 
tion of the interface and on the type of scattering, specular or 
diffusive, which determines whether an electronic trajectory 
crosses both bands or only one. Thus, knowledge of the junc- 
tion geometry and interface properties should make it possible 
to compare the model with experimental data. Qualitatively, 
one can see already that the observation of a zero-bias conduc- 
tance peak can be consistent with nodeless superconductivity, 
and that a non-zero energy surface bound state, that can exist 
at subgap as well as supergap energies, can easily be mistaken 
for a gap feature when interpreting conductance spectra. 

This work was partially supported by the Netherlands Or- 
ganization for Scientific Research (NWO) and the NanoNed 
program under Grant No. TCS7029. 



FIG. 3: Bound state energy at contacts to s± superconductors as 
function of the band ratio parameter a for A2 = 2Ai. No bound 
states exist at ene rgies bet ween the two gaps (shaded region). Also, 
for band ratios y^AT/A^ < a < 1 no surface bound states exist. 



5 



4 

3 

> 









1 1 1 1 1 

Z=10, A 2 =2A 1 








a=0 




a=0.2 




a=0.5 




a=1 




I a=10 












_>^l^*v 












/' ; 







12 3 4 

eV/A, 



FIG. 4: (Color online) Conductance in the low transparency regime, 
Z — 10, in the s++ model. 

Still, interference effects at a ~ 1 result in a complex dl/dV 
behavior, where the conductance is not a simple sum over two 
individual bands. Presently, the conductance spectra of con- 
tacts with the multiband superconductor MgB 2 are usuall y fit - 
ted by the sum of two single-band tunneling probabilities II 1 411 . 
With the increased level of sophistication in the realization of 
epitaxial magnesium diboride junctions and single crystalline 
point contacts, one can expect that the present predictions for 
multiband interference effects can be observed there as well. 
In order to demonstrate the main features of the model, 
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